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We show that bright sohtons exist in quasi-one dimensional heteronuclear multicomponent Bose- 
Einstein condensates with repulsive self-interaction and attractive inter-species interaction. They 
are remarkably robust to perturbations of initial data and collisions and can be generated by the 
mechanism of modulational instability. Some possibilities for control and the behavior of the system 
in three dimensions are also discussed. 
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I. INTRODUCTION 

Symbiosis is an assemblage of distinct organisms living 
together. Although the original definition of symbiosis 
by De Bary 1] did not include a judgment on whether 
the partners benefit or harm each other, currently, most 
people use the term symbiosis to describe interactions 
from which both partners benefit. 

In Physics, waves in dispersive linear media tend to 
expand due to the different velocities at which the wave 
components propagate. This is not the case in many non- 
linear media, in which certain wavepackets, called soli- 
tons are able to propagate undistorted due to the balance 
between dispersion and nonlinearity 0. 

Stable solitons of different subsystems are sometimes 
able to "live together" and form stable complexes called 
vector solitons as it happens with Manakov optical soli- 
tons or stabilized vector solitons 5]. In some cases, 
a (large) robust soliton can be used to stabilize a (small) 
weakly unstable wave 6]. 

Multicomponent solitary waves also appear in Bose- 
Einstein condensates (BECs). In fact, multicomponent 
BECs support nonlinear waves which do not exist in sin- 
gle component BECs such as domain wall solitons 
dark-bright solitons 9], etc. Most of the previous analy- 
ses correspond to homonuclear multicomponent conden- 
sates for which the atom-atom interactions are repulsive. 
However, heteronuclear condensates offer a wider range 
of possibilities, the main one being the possibility of hav- 
ing a negative inter-species scattering length. This pos- 
sibility has been theoretically explored in the context of 
Feschbach resonance management IQjJ and realized ex- 
perimentally for boson- fermion mixtures 0, . 

In this paper we study the existence and properties 
of bright solitons in heteronuclear two-component BECs 
with scattering lengths aii,a22 > and ai2 < 0. We 
would like to stress the fact that these coefficient combi- 
nations do not arise in other systems where similar model 
equations are used. For instance in nonlinear optics, 
where the nonlinear Schrodinger equations used to de- 
scribe the propagation of laser beams in nonlinear media 
are similar to the mean field equations used to describe 
Bose-Einstein condensates, the nonlinear coefficients are 
allways of the same sign. The closest analogy could 



happen in the so-called QPM ( quasi-phase-matched ) 
quadratically nonlinear media, where an effective cubic 
nonlinearity could be "engineered" which could have sim- 
ilar properties but we do not know of any systematic 
studies of those systems. 

Our analysis will show novel features with respect to 
those already found in single species BECs lm. For 
instance, even when solitons do not exist for each of 
the species, the coupling leads to robust vector solitons. 
Since the mutual cooperation between these structures 
is essential for their existence we will refer to these soli- 
tons hereafter as symbiotic solitons. We also show how 
they appear by modulational instability and study some 
features of their collisions. We also comment on the pos- 
sibility of obtaining these structures in multidimensional 
configurations. 



II. THE MODEL AND ITS BASIC PROPERTIES 

In this paper we will study two-component BECs in 
the limit of strong transverse confinement ruled by 
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where x is the adimensional longitudinal spatial variable 
measured in units of ao = \/h/miUJ^^ t is the time mea- 
sured in terms of l/cc;^, and Uj{xA) = Uj{r^ r) 
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with a = mi/m2 and a^j being the 5— wave scattering 
lengths. The normalization for Uj is j\uj\^d^x = Nj 
where Nj is the number of particles of each species. 

Let us first consider constant amplitude solutions of 
Eq. (QJ, which are of the form 
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for j = 1,2. We will study the evolution of small pertur- 



bations of (pj of the form 
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Using Eq. (QJ and retaining the first order terms we get 
partial differential equations for SAi , S/3i , SA2 , S/32 which 
can be transformed to Fourier space to obtain 



(4) 



ao{k) being the Fourier transform of the initial pertur- 
bation. Perturbations remain bounded if Re [ri(fc)] < 0. 
Some algebra leads to 



= i (/l +/2 ± v/(/l-/2)2 + 4C2) 
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where f, = -{g,,Aj + e/A)k\C^ = AfA^gf^k^ 
The so-called modulational instability (MI) occurs when 
ft{k)'^ > for any k. For small wavenumbers (worst sit- 
uation) we get 



9l2 > 911922, 
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which is analogous to the miscibility criterion for two- 
component condensates Q. However, the physical mean- 
ing of Eq. is very different since now this instability 
is a signature of the tendency to form coupled objects 
between both atomic species. The role of MI in the for- 
mation of soliton trains and domains in BEC has been 
recognized in previous papers 0, 0, • 



III. VECTOR SOLITONS 



Eqs. (pP) have sech-type solutions 
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and the MI condition are satisfied. Eq. implies 
that, given the number of particles in one component the 
other is fixed. 

Since the self-interaction coefficients are positive, these 
solitons are supported only by the mutual attractive in- 
teraction between both components. This type of vector 
soliton thus differs from others described for Nonlinear 
Schrodinger equations of the form Eq. such as the 
Manakov solitons "s] , where all the nonlinear coefficients 
cooperate to form the solitonic solution. The MI condi- 
tion implies that the formation of these solitons has 
a threshold in gi2 and means that the cross-interaction 
must be strong enough to be able to overcome the self- 
repulsion of each atomic cloud. There are no analogues 
to this condition in single component systems since soli- 
tons exist for any value of the self-interaction coefficient 
^ < 0. To fix ideas, taking a ^^Rb-^^K mixture with 
ail = 69ao,a22 = 99ao the MI condition implies that 
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FIG. 1: Dependence of the ratio N2/N1 of sech-type vector 
solitons on the inter-species scattering length ai2. (a) 
^^K mixture with an = 69ao,a22 = 99ao, (b) An hypotetical 
^Li-^^Na mixture with an = 5ao and a22 = 52ao. 
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FIG. 2: [Color online] Evolution of displaced soliton initial 
data of the form ui = {Ni/{2w)y^'^sech{{x -\- xo)/w), U2 = 
{N2/{2w)y^'^sech(x/w) for a ^^Rb-^^K mixture with Ni = 
3000, ATs = 1189, an = 69ao,ai2 = -90ao,a22 = 99ao in a 
trap with (jtj± — 215 Hz. 



ai2 < — 83ao in order to obtain solitons. In Fig. fflja) it 
can be seen how the ratio N2/N1 is close to 0.4 in the 
range of values of — 83ao > ai2 > — ISOao- An hypothet- 
ical ^Li-^^Na mixture with an = 5ao and a22 = 52ao (in 
appropriate quantum states) leads to the curve in Fig. 
^b), which shows a much larger range of variation. 



IV. SOLITON STABILITY 

We can use the Vakhitov-Kolokov (VK) criterion to 
study the stability of solitons given by Eq. {TJ. To do 
this, we must study the sign of dXj/dNj. For soliton 
solutions this can be done from the explicit form of Xj. 
After some algebra we find Xi{Ni) and X2{N2) and ob- 
tain that dXi/dNi > 0, and SA2/S7V2 > in ah their 
range of existence, which proves the linear stability of 
the solitons for small perturbations and contradicts the 
naive intuition that the self-repulsion would lead to in- 
trinsically unstable wavepackets. 

We have studied numerically the robustness of sym- 
biotic solitons to finite amplitude perturbations. First 
we have perturbed both solutions with small amplitude 
noise and found that, in agreement with the predictions 
of the VK criterion, they survive after the emission of 
the noise in the form of radiation. Next we have applied 
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FIG. 3: [Color online] Evolution of the ground state of a 
^'^Rb-'^^K mixture with iVi = 25000, iV2 = 20000 after switch- 
ing the interspecies scattering length from ai2 = 95ao to 



ai2 = — 90ao. (a) Initial state: 



(blue dotted line) and 



1 2/2 1 (red solid line), (b) Profile of for t = 500 show- 

ing three remaining solitons. (c) Pseudocolor plot of for 
X e [-500, 500] and t G [0, 500]. 



a stronger perturbation consisting of displacing mutually 
their centers and observe that a soliton is formed even 
for relative displacements of the order of the soliton size 
[Fig. 121. Finally we have started with sech-type initial 
data which are not solitons and observe that after the 
emission of some radiation solitons are formed. 



V. GENERATION OF SYMBIOTIC SOLITONS 
BY MI 

To study the generation of these solitons by MI in real- 
istic systems we have considered a multicomponent Bose- 
Einstein condensate of ^^Rb and ^^K atoms for which the 
inter-species scattering length ai2 is controlled by the use 
of Feschbach resonances as proposed in . To simplify 
the problem here we do not consider the effect of gravity. 

We start by constructing the ground state of the sys- 
tem for an elongated trap typical of the LENS setup 
with = 215 Hz, CO = 16.3 Hz. For these atomic 
species an = 69ao and a22 = 99ao. We adjust the inter- 
species scattering length to ai2 = 95ao during the con- 
densation process. The ground state of this system for 
A^i = 25000, N2 = 20000, shown in Fig. EJa), agrees well 
with the theoretical predictions for these systems ^2^. 

After the condensate is formed we change instanta- 
neously this quantity to a negative value and at the same 
time switch off the longitudinal trapping potential and 
observe numerically the evolution of the ground state. 

First we choose ai2 = — 90ao and observe the evolution 
starting from the ground state with ai2 = 95. Since 
the inter-component repulsive force is not present now, 
the sharp domain wall separating both species (see Fig. 
OJa)) decay through a highly oscillatory process related 
to the formation of a shock wave 0| . The final outcome 
is the formation of a soliton train (see Fig. |3Jb,c)) of 
which three solitons of about 20 jam size and each with 
about 3000 rubidium and 1200 potassium atoms remain 
in our simulation domain after 500 adimensional time 




FIG. 4: [Color online] Pseudocolor plots of the evolution of the 
ground state of a ^'^Rb-^^K mixture with iVi = 25000, A^2 = 
20000 after switching the inter-species scattering length from 
ai2 = 95ao to: (a) ai2 = — 70ao (below the MI limit); (b) 
ai2 = — 87ao (slightly above the MI limit a* = — 83ao. 



units [Fig. ISI^c)]. Other smaller and wider solitons exit 
our integration region traveling at a faster speed. 

The final number of solitons depends on the value of 
ai2 choosen during the condensation process (which con- 
trols the overlapping of the species) and the number of 
particles Ni^ N2 and the negative scattering length ai2 
choosen to destabilize the system. For instance, choos- 
ing ai2 = — 70ao, which is below the theoretical limit 
for MI the evolution of the wavepacket is purely disper- 
sive [see Fig. ^a)]. Choosing ai2 = — 87ao, above the 
MI limit but below the choice of Fig. |31 leads to the 
formation of a single soliton [Fig. E^b)]. It seems that 
the larger the scattering length, the larger the number of 
solitons which arise after the decay of the initial config- 
uration. The many degrees of freedom present in these 
system open many posibilities for controling the number 
and sizes of solitons by appropriately choosing the values 
of ai2 before and after the condensate is released and the 
initial number of particles Ni,N2. 



VI. COLLISIONS OF SYMBIOTIC SOLITONS 

The robustness of symbiotic solitons manifests also 
in their collisional behavior and their internal structure 
makes the interaction of these vector solitons very rich. 
Since each soliton is a compound object the collisions are 
at the same time a coherent phenomenon because of the 
direct overlapping of the same type of atoms and an in- 
coherent one because of the incoherent nature of interac- 
tion between different types of atoms. A related subject 
of recent interest in Optics is that of partially coherent 
solitons 23]. 

We have simulated head-on collisions of equal symbi- 
otic solitons of opposite velocities given by 




for j = 1,2. aj^± are the relative phases and N2 is 
given by Eq. In Fig. [SI we show some examples 
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FIG. 5: [Color online] Head-on collisions of symbiotic solitons 
with A^i 3000, A^2 1189, 'u; 1.723, an 69ao,ai2 
— 90ao and a22 = 99ao. (a)-(d) Slow collisions for v — 0.05 
and (a-b) = 0^2,+ , 0^2,-) = (0,0,0,0), (c) 6 = 

(0,7r,0,0), (d) = (7r/2, 0,0,0). Moderate speed collisions 
(e) V = 0.2, 6 = (0, 0, 0, 0), (f) v = 0.2, 6 = (7r/2, 0, 0, 0). 



of these collisions. Slow [Fig. [Sfa-d)] or moderate speed 
collisions [Fig. ISte-f)] lead to bound solitons while for 
larger speeds the picture is not so clear. The specific 
outcome of the collision depends on the relative soliton 
phases with the phase difference between the larger com- 
ponents in the symbiotic soliton (in this case Rb) being 
the dominant ones. For instance collisions with phases 
= (ai,+, a2,+, ^2,-) = (0,7r,7r,0) [Fig. ISl(c)] and 
= (0, 0, TT, 0) (not shown) both lead to mutual repulsion 
but the outgoing speeds are different due to the differ- 
ent interactions between the internal components of the 
soliton. Collisions with higher but still moderate speeds 
[Fig. inije-f)] give independent vector solitons. The out- 
come of the collisions with zero phases is a bound state 
of two vector solitons which has internal oscillations, i.e. 
some sort of multicomponent higher order soliton. 



VII. PROSPECTS FOR MULTIDIMENSIONAL 
SYMBIOTIC SOLITONS 

A very interesting question arising naturally is: do 
these symbiotic solitons exist in multidimensional sce- 



narios? In principle the answer is not evident since the 
only effect acting against stabilization of multidimen- 
sional soliton structures would be collapse, but one could 
think that in this case collapse could be inhibited because 
of the repulsive self-interaction, thus a deeper analysis is 
in order. 

The adimensional model equations in two and three 
dimensions take the form 
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with j = 1, 2 and /c = 2, 1 correspondingly. 

Let us first consider this problem in two spatial dimen- 
sions. To study collapse rigorously one usually tries to 
compute the exact evolution of the wavepacket widths 
rigorously 24]. For the multicomponent case and mi = 
m2 = m, this was studied by group-theoretical methods 
by In our case, from the general formulae obtained 
by Gosh we get a sufficient condition for collapse, which 
is 
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In principle, this is a bad result for obtaining localized 
structures since it means that arbitrarily close to any 
stationary solution (for which = 0), there would be 
collapsing solutions and thus stationary solutions, if they 
exist, would be unstable. As it is usual in the framework 
of collapse problems the situation would be even worse 
in three spatial dimensions with solutions of arbitrary 
small number of particles undergoing collapse provided 
they are initially sufficiently localized. 

This means that in principle symbiotic solitons could 
only be obtained in quasi- ID geometries because of the 
transverse stabilization effect provided by the trap in a 
similar way as ordinary bright solitons do. 



VIII. CONCLUSIONS AND EXTENSIONS 

In this paper we have studied vector solitons in het- 
eronuclear two-component BECs which are supported by 
their attractive mutual interaction. These symbiotic soli- 
tons are linearly stable and remarkably robust and can be 
generated through modulational instability phenomenon 
with many possibilities for control. Collisions of these 
vector solitons show their robustness and open different 
ways for their manipulation and the design of novel quan- 
tum states such as breather-like states. We have also con- 
sidered multidimensional configurations and shown that 
collapse may avoid the formation of fully multidimen- 
sional symbiotic solitons. 

We think that the conceptual ideas behind our work 
can also be used to understand boson-fermion mixtures. 
For instance, ai2 is known to be negative and large for 
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quantum degenerate mixtures of ^''Rb and 0. In 
those systems numerical simulations have proven the for- 
mation of localized wavepackets which could share 
the same essential mechanisms for the formation of soli- 
tary waves. 
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